Collinear to Anti-collinear Quantum Phase Transition by Vacancies 
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We study static vacancies in the collinear magnetic phase of a frustrated Heisenberg J1-J2 model. 
It is found that vacancies can rapidly suppress the collinear antiferromagnetic state (CAFM) and 
generate a new magnetic phase, an anti-collinear magnetic phase (A-CAFM), due to magnetic 
frustration. We investigate the quantum phase transition between these two states by studying a 
variety of vacancy superlattices. We argue that the anti-collinear magnetic phase can exist in iron- 
based superconductors in the absence of any preceding structural transitions and an observation 
of this novel phase will unambiguously resolve the relation between the magnetic and structural 
transitions in these materials. 
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There are several reasons for studying static vacancy 
problems on frustrated magnetic systems. First of all, 
there has been convincing experimental evidence which 
supports the magnetism in iron-based high temperature 
superconductors (Fe-HTSC) can be understood by an 
effective frustrated magnetic model {Ji-J%-Jz model) [H- 
IH which simultaneously captures the collinear antiferro- 
magnetic state and the tetragonal to orthohombic struc- 
tural transition observed in neutron-scattering experi- 
ments Q. The new superconductors are very flexible 
in substituting Fe by other transition metal atoms, such 
as Mn, Zn, Co and Ni. The static-vacancy problem in 
the 3\-Ji-J z model is, then, an important low energy ef- 
fective model for non-magnetic Zn-doped Fe-HTSC Q. 
Moreover, the recently discovered 122 iron-chalcogenide, 
(K,Cs)Fe2- x Se2 |2H9|, carries intrinsic iron vacancies, 
which can even form superlattice vacancy structures [lol - 
[Hj]. Thus, the solution of the static- vacancy problem 
can be directly tested experimentally and contributes to 
a fundamental understanding regarding the role of mag- 
netism in superconductivity as well as the coupling be- 
tween lattice and magnetism. Second, with various frus- 
trated magnetic materials being discovered in the past 
decade, many novel physics and new states of matter 
have been proposed. However, experimentally, it has of- 
ten been difficult to identify features associated to novel 
physics, for example, spin liquid state [16|, ll7(. Static va- 
cancies can either enhance or decrease the degree of frus- 
tration and can behave rather differently in different state 
of matters. Therefore, static vacancies can contribute to 
a new understanding of frustrated magnetic physics and 
provide unique features that can be probed experimen- 
tally. Finally, even in a standard quantum Heisenberg an- 
tiferromagnetic model, it has been shown that quantum 
fluctuations can also be dramatically modified around 
static vacancies [3, Studying static vacancies in 

frustrated quantum magnetic systems can also provide a 
deeper understanding of the interplay between quantum 
fluctuations and geometric frustration. 



In this Letter, we study the static vacancy problem in 
the J1-J2 antiferromagntic Heisenberg model. We em- 
ploy a linear spin-wave (LSW) theory 18] to understand 
properties of a single static vacancy and static vacancy 
superlattices. We show, depending on the frustrated cou- 
pling, quantum fluctuations can be either reduced or en- 
hanced on neighbors of an isolated vacancy. More im- 
portantly, by calculating the exact ground-state prop- 
erties of a variety of static vacancy lattices, we predict 
that sufficient static vacancies can cause a quantum phase 
transition between the collinear magnetic phase and an 
anti-collinear magnetic phase before a spin glassy phase 
without a spatial long-range magnetic order is formed. 

Without vacancies, the J\-Ji model is given by 



Hq — J\ Si ■ Sj + J2 Si ■ Sj, 



(1) 
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where < ij >nn and < ij >nnn denote bonds formed 
by two nearest neighbor sites and two next nearest neigh- 
bor sites respectively. For a classical J1-J2 model with 
Ji < 2J2, the ground state can be viewed as two decou- 
pled antiferromagnetically (AFM) ordered states on the 
A and B sublattices as shown in Fig[TJ Including quan- 
tum fluctuations, the relative angle between the two anti- 
ferromagnetic orders on the A and B sublattices is locked 
and the quantum model has a CAFM ground state with 
an ordering wave vector at Q = (0, n) or Q' = (tt, 0). The 
CAFM state is driven by the frustrated coupling J\ . The 
energy of quantum fluctuations can be calculated using 
the standard LSW theory. Without losing generality, we 
take the AFM order in the A sublattice as S Z A ^ and 
the AFM order in the B sublattice rotates by 9 around 
y-axis relative to the one in the A sublattice. Namely, 
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In the LSW approximation, Eq.(p} reduces to 
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H = Y,"kb{b k + f(b_ k b k + h.c.)-4NJ 2 S 2 , (3) 
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FIG. 1: (color online) The sketch of the CAFM state in 
the Ji — J2 model where the two sublattices A and B are 
colored by black and red respectively and a single vacancy 
site is at the center. The numbers labels the magnetic order 
parameter < S z (i) > at each site in the CAFM state in the 
presence of a single vacancy for the parameters J\ = J2 and 
5=1 (without the vacancy, the magnetic order parameter is 
I < S z (i) > I = 0.7817). 



where, a;*; = S[Jicos8(cosk x —cosk y )+Ji(cosk x +cosk y )+ 
4J2], v k = S[JiCos6(cosk x — cosky) — ^(coskx + cosky) — 
4J2COsfc x cosfc„], and b (b + ) are magnon annihilation (cre- 
ation) operators. Using the Bogliubov transformation, 
b k = cosh-0feo;fe + sinh-0feQ;l fc , Eq.([3]) can be diagonalized 
as 



H = Y;Uko{a k + Ne - 4NJ 2 S{S + 1), 



(4) 



where, ui k = y/u% - v\, 

cosh 2 ?/; = | 
and 
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, smh 2 ?/; = i 



- 1 



N 9 



./ ojI — vl = c — an 
N ^ 2 V k J 2 

k 



r2o 

^COS 2 0, (5) 



where c is independent of 9 and clq ~ 0.033. eo describes 
the well-known "order by disorder mechanism" and has a 
minimum at 9 = 0, tt, favoring a CAFM order [l, 2(J 21 1. 

For a single vacancy at the origin of the A-sublattice, 
the total Hamiltonian can be written as H = Ho — V 
with V = Ji S ■ Sf + J 2 V. So ■ * A 

LSW approximation, the Hamiltonian becomes 



Sf. In the 
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k+q^k 



-k+qCt k ) 
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where A fc , 9 = Ak^cosh^k+q + V>fc) + z^sinh^/c+g + 
ipk), B k . q = A kt qsmh.(ip k +q + ipk) + v k cosh(ip k+q + 
-0fe), an d C k = p k (coship k + sidhipk) with p k = 
—2SJ\ S/2sin9(cosk x — cosk y ). Up to the first order 
of J 2 /J2, the total ground state energy in the presence 
of a single vacancy is E = -4(iV - 1) J 2 S(S + 1) + (N - 
l)eo + e v , where 
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(7) 



where a\ ~ 0.36. e u has a minimum at = ±7r/2 and 
does not favor a CAFM state. The physics behind the 
energy e v can be argued as follows. In the CAFM state, 
creating a vacancy at one sublattice is similar to applying 
an external magnetic field along magnetic ordered direc- 
tion on the four neighbor sites of the vacancy in the other 
sublattice. Since the spins of the four neighbor sites are 
AFM, the presence of such a field would favor the AFM 
order in the four neighbor sites to be perpendicular to the 



external magnetic field direction. e and e v have differ- 
ent dependence on the spin S. The competition between 
these two energies can lead to a new phase transition. 

Considering the model with a small density of vacan- 
cies, p, in the first order approximation and up to a 
constant, we can approximate the energy density of the 
model as a function of 9 to be 



e(9,p) = (l-p)eo(9)+pe v (9). 



(8) 



The energy density favors the CAFM state(0 = 0, tt) if 
p < p c and an A-CAFM state (9 = ±7r/2) if p > p c where 
the critical vacancy density is given by 



Pc = 



do 



a + a\S 



(9) 



Pluging in the values of clq and 01, we obtain p c — 0.086 
for S = 1 and p c = 0.158 for 5=1/2. These criti- 
cal values are well below the percolation threshold which 
destroys the long range AFM order. 
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We can also solve the single vacancy problem exactly 
(within the LSW approximation). Defining the standard 
Green functions: 



G j , j ,(t) = -i<T[b j (t)b+(0)] > 
F jd ,(t) = -i<T[bf(t)b+(p)] >, 



(10) 



can 



and their Fourier transformation G(F)j,j> (t) 

E q e^e-^'-^ + ^G(F) w , we 
derive the following dynamic equations for the Green 
functions in the presence of a single vacancy at the 
origin of the lattice, 
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Jj [Ak+q,p{G° k+q G p , k + F k+g F Pik ) + B k+qtV (G° k+q F p>k + F k+q G Pik )] , 
p 

"jy 51 [Ak+g^iGl+qFp^ + F" +q G p .k) + B k+q ^(G k+q Gp^ k + F k+q F Pyk )] , 



(11) 



where G(F)° are given by 
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and A k + q , p = J X S . 

■ k„ — p y )] + cos#[cos(g : 



(12) 
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■ 4J 2 5[1 + cos(g : 
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2cos6*[cos(q. E + k x 
)] + cos#[cos((fc + k x ) — cos(q y + k y ) - 
[cos(q x + k x ) + cos(q y + k y ) + cosp x - 

+ cos(q x + k x - p x )cos(q y + k y - p y )\ 



B k 



J\S\ cos0[cos{q x + kx) — cos(q y + k y ) + COSpa; 



k+q,p 

cosp y ] — [cos(q x + k x ) + cos(q y + k y ) + cosp x + cosp a ]J — 

4 J 2 S'[cosp :E cosp a + cos(q x + k x )cos(q y + k y )]. The Dyson- 
type equations in Eqs. tjlip can be solved numerically for 
any given 9 and J\jJ% values. We focus on the mag- 
netic order moments and the total energy on sublattices 
surrounding the vacancy located at the origin (0,0). 

First, in Figfl] we report the magnetic order parame- 
ter < S z (i) > at each site in the CAFM state (9 = 0) for 
the parameters J\ — J 2 and 5 = 1 (without the vacancy, 
the uniform magnetic order is | < S z {i) > | = 0.7817). In 
FigfTJ we plot the magnetic moments at the sites (0, 1), 
(1, 0) and (1, 1) as a function of Jx/J-z in the CAFM state. 
There are two important results: (1) the effects of the 
vacancy on its nearest neighbor (NN) sites are different 
along the two directions in the CAFM state. The zero- 
point deviations are suppressed (enhanced) at the NN 
sites along the ferromagnentic (AFM) directions if J\ is 
AFM and the results reverse if J\ is negative (ferromag- 
netic); (2) the effect of the vacancy on its next nearest 
neighbor (NNN) does not break C4 rotation symmetry 
even in the CAFM state. The zero-point deviation at 
these sites is suppressed for small |Ji| values. This re- 
sult is not surprising since it is known to be the case for 
J\ = 0. However, the deviation goes from depression 
to enhancement as | Ji| increases further. This crossover 
reflects the frustration increases the transverse fluctua- 
tions due to the anti-collinear tendency between the two 
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FIG. 2: (color online) The ^-dependence of energy per spin 
in the presence of single vacancy for three different size of 
clusters: 3 x 3, 5 x 5 and 7x7. The parameters are chosen 
as J2 — Ji — 1 and S = 1. 



magnetic moments of the sublattices around the vacancy. 

Second, we calculate the total energy of the model on 
clusters centered at the static vacancy as a function of 
6. In Fig(21 we plot the energy on three different clus- 
ters surrounding the vacancy with sizes, 3x3, and 5x5 
and 7x7 and parameters S — 1, J\ — Ji. It is clear 
that the energy minimum for a 3 x 3 cluster is 9 = tt/2. 
Moreover, the static magnetization at the 8 sites in the 
3x3 lattice is around 0.9/i& for 9 = tt/2 which is larger 
than the case with no vacancies 0.798/is at the CAFM 
phase. This result confirms that the vacancy clearly fa- 
vors an A-CAFM ordering between two sublattices. In 
Figj3ja), we plot the configuration of magnetic moment 
surrounding the vacancy in the A-CAFM state 9 = tt/2 
with Jx = J2 = 1 and 5 = 1. The value of the magnetic 
moment along the z direction for the nearest neighbour 
site of the vacancy linearly increases as a function of J\ 
as shown in Figj3jb). 

The above study of a single vacancy suggests that the 
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FIG. 3: (color online) (a) The magnetic moment configura- 
tion around the vacancy in 3 x 3 cluster with Ji = Ji = 1 and 
5=1. (b) The Ji-dependence of ) and S^ ^ for S = 1 
and J2 = 1. 




FIG. 4: (color online) The energy of three different superlat- 
tices as a function of 8 for J\ = J%. For the 4x4 lattice, the 
energy minimum becomes 9 — 0. 



A-CAFM configuration is favored if only the energy on 
the small sublattice surrounding the vacancy is consid- 
ered. In order to confirm that the existence of the global 
phase transition in the presence of vacancies, we take a 
super unit cell in the square lattice with N x M sites 
and creates one vacancy in the unit. Thus, if we repeat 
this unit to create a superlattice, we obtain a system in 
which the percentage of vacancy concentration is given by 
N 1 M ■ In this superlattice system, for a given wavevec- 
tor k, the Eq. ([TTj) can be reduced to equations that only 
couple 2N x M Green functions given by Gk+Q nm ,k 
and Fk+Q nm ,k, where Q n ,m = (2irn/N,2irm/M) and 
n(m) = 0, N(M) — 1. In FigJH we show the energy 
of three different supperlattices as a function of 8 for 
Ji = J2 and 5 = 1. For both superlattices with 2x4 
and 2x6 unit cells which are corresponding to 12.5% 
and 8.4% vacancy concentration respectively, the anti- 
collinear state is favored. However, the CAFM state is 
favored in a superlattice with 4x4 unit cell corresponding 
to 6.3% vacancy concentration. This result justifies our 
previous rough estimation of the critical vacancy density. 

Our above results have important implications in iron- 
based superconductors. All of our above calculations 



demonstrate an existence of quantum phase transition 
from a CAFM state to an A-CAFM state at a certain 
critical vacancy concentration p c . While the CAFM state 
breaks C 4 rotation symmetry, the A-CAFM state does 
not break C 4 rotation symmetry. In iron-pnictides, there 
is always a tetragonal-to-orthorhombic structural transi- 
tion which occurs at the temperature above or equal to 
CAFM transition temperature. This structural transi- 
tion breaks C 4 to C 2 and is naturally explained as a 
consequence of magnetic fluctuations associated with the 
CAFM state Q, H[. If the A-CAFM state exists and the 
structural transition is magnetically driven, our results 
predict that the lattice distortion can be absent in the 
A-CAFM phase. 

It is also worth to discuss that the vacancy or- 
derings have been observed in (Ax- y Fe2- x Se2) iron- 
chalcogenides, where the vacancy patterns are corre- 
sponding to a natural reduction of the magnetic frus- 
tration so that the magnetic transition temperature is 
strongly enhanced [22j. The vacancy superlattices used 
in our calculation do not reduce the magnetic frustra- 
tion. Therefore, our results do not directly apply to the 
observed vacancy patterns, such as the 245 pattern in 
K 2 FeiSe5 flol | . However, for the materials with very di- 
luted vacancy concentration, we expect that our result 
should be valid as well. 

In summary, we study static vacancies in the collinear 
magnetic phase of a frustrated Heisenberg J1-J2 model 
and identify a quantum phase transition between 
collinear antiferromagnetic state (CAFM) and an anti- 
collinear antiferromagnetic phase (A-CAFM). Our re- 
sults can help to resolve the relation between magnetic 
and structural transitions in iron-based superconductors. 
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